Uncovering Protein Ensembles: Automated Multiconformer Model Building for X-ray Crystallography and Cryo-EM.

In their folded state, biomolecules exchange between multiple conformational states that are crucial for their function. Traditional structural biology methods, such as X-ray crystallography and cryogenic electron microscopy (cryo-EM), produce density maps that are ensemble averages, reflecting molecules in various conformations. Yet, most models derived from these maps explicitly represent only a single conformation, overlooking the complexity of biomolecular structures. To accurately reflect the diversity of biomolecular forms, there is a pressing need to shift towards modeling structural ensembles that mirror the experimental data. However, the challenge of distinguishing signal from noise complicates manual efforts to create these models. In response, we introduce the latest enhancements to qFit, an automated computational strategy designed to incorporate protein conformational heterogeneity into models built into density maps. These algorithmic improvements in qFit are substantiated by superior Rfree and geometry metrics across a wide range of proteins. Importantly, unlike more complex multicopy ensemble models, the multiconformer models produced by qFit can be manually modified in most major model building software (e.g. Coot) and fit can be further improved by refinement using standard pipelines (e.g. Phenix, Refmac, Buster). By reducing the barrier of creating multiconformer models, qFit can foster the development of new hypotheses about the relationship between macromolecular conformational dynamics and function.

Refinement of multiconformer ensemble models from multi-temperature X-ray diffraction data.

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363 K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Refinement of Multiconformer Ensemble Models from Multi-temperature X-ray Diffraction Data.

Conformational ensembles underlie all protein functions. Thus, acquiring atomic-level ensemble models that accurately represent conformational heterogeneity is vital to deepen our understanding of how proteins work. Modeling ensemble information from X-ray diffraction data has been challenging, as traditional cryo-crystallography restricts conformational variability while minimizing radiation damage. Recent advances have enabled the collection of high quality diffraction data at ambient temperatures, revealing innate conformational heterogeneity and temperature-driven changes. Here, we used diffraction datasets for Proteinase K collected at temperatures ranging from 313 to 363K to provide a tutorial for the refinement of multiconformer ensemble models. Integrating automated sampling and refinement tools with manual adjustments, we obtained multiconformer models that describe alternative backbone and sidechain conformations, their relative occupancies, and interconnections between conformers. Our models revealed extensive and diverse conformational changes across temperature, including increased bound peptide ligand occupancies, different Ca2+ binding site configurations and altered rotameric distributions. These insights emphasize the value and need for multiconformer model refinement to extract ensemble information from diffraction data and to understand ensemble-function relationships.

Ligand binding remodels protein side-chain conformational heterogeneity.

While protein conformational heterogeneity plays an important role in many aspects of biological function, including ligand binding, its impact has been difficult to quantify. Macromolecular X-ray diffraction is commonly interpreted with a static structure, but it can provide information on both the anharmonic and harmonic contributions to conformational heterogeneity. Here, through multiconformer modeling of time- and space-averaged electron density, we measure conformational heterogeneity of 743 stringently matched pairs of crystallographic datasets that reflect unbound/apo and ligand-bound/holo states. When comparing the conformational heterogeneity of side chains, we observe that when binding site residues become more rigid upon ligand binding, distant residues tend to become more flexible, especially in non-solvent-exposed regions. Among ligand properties, we observe increased protein flexibility as the number of hydrogen bonds decreases and relative hydrophobicity increases. Across a series of 13 inhibitor-bound structures of CDK2, we find that conformational heterogeneity is correlated with inhibitor features and identify how conformational changes propagate differences in conformational heterogeneity away from the binding site. Collectively, our findings agree with models emerging from nuclear magnetic resonance studies suggesting that residual side-chain entropy can modulate affinity and point to the need to integrate both static conformational changes and conformational heterogeneity in models of ligand binding.

Proteins are the workhorses of our cells. They are large molecules that 'fold' into specific, often highly complex, three-dimensional configurations. These structures are not static, but rather dynamic and flexible. In other words, proteins can shift between different three-dimensional shapes to perform their tasks within the cell. To perform their roles, many proteins have to bind to small molecule ligands. Many ligands are drugs, which means that their effectiveness depends on their ability to bind to and impact the proteins involved in the disease they are treating. When a ligand binds to a protein, it can reshape the protein. For example, certain conformations of the protein, which were difficult for the protein to be in on its own, may become more stable when the ligand binds. Additionally, upon ligand binding, some parts of the protein may move relative to each other. Previous studies have shown that these movements can affect the interaction between ligand and protein. However, these studies only examined a small number of proteins. Therefore, Wankowicz et al. set out to determine, in greater detail, what happens to protein flexibility upon ligand binding. First, a pipeline was created to model alternative configurations of the protein both with and without ligands attached. These models measured flexibility within protein structures. The models revealed that when ligands bind to proteins, the flexibility of different regions of the protein changes ­ and does so in a consistent way. Proteins that become more rigid in the region interacting with their ligands become less rigid in other, distant regions, and vice versa. In other words, the rest of the protein is able to compensate for any changes in flexibility caused by ligand binding, which may contribute to how well a ligand binds to a protein. This study demonstrates the ability of ligands to affect the entire structure of the proteins they bind to, and therefore sheds new light on the role of proteins' innate conformational flexibility during this process. These results will contribute to our understanding of how the ligands and proteins involved in different cellular processes interact with each other ­ and, potentially, how these interactions can be manipulated.

Molecular features of exceptional response to neoadjuvant anti-androgen therapy in high-risk localized prostate cancer.

High-risk localized prostate cancer (HRLPC) is associated with a substantial risk of recurrence and disease mortality. Recent clinical trials have shown that intensifying anti-androgen therapies administered before prostatectomy can induce pathologic complete responses or minimal residual disease, called exceptional response, although the molecular determinants of these clinical outcomes are largely unknown. Here, we perform whole-exome and transcriptome sequencing on pre-treatment multi-regional tumor biopsies from exceptional responders (ERs) and non-responders (NRs, pathologic T3 or lymph node-positive disease) to intensive neoadjuvant anti-androgen therapies. Clonal SPOP mutation and SPOPL copy-number loss are exclusively observed in ERs, while clonal TP53 mutation and PTEN copy-number loss are exclusively observed in NRs. Transcriptional programs involving androgen signaling and TGF-ß signaling are enriched in ERs and NRs, respectively. These findings may guide prospective validation studies of these molecular features in large HRLPC clinical cohorts treated with neoadjuvant anti-androgens to improve patient stratification.

Cryo-EM model validation recommendations based on outcomes of the 2019 EMDataResource challenge.

This paper describes outcomes of the 2019 Cryo-EM Model Challenge. The goals were to (1) assess the quality of models that can be produced from cryogenic electron microscopy (cryo-EM) maps using current modeling software, (2) evaluate reproducibility of modeling results from different software developers and users and (3) compare performance of current metrics used for model evaluation, particularly Fit-to-Map metrics, with focus on near-atomic resolution. Our findings demonstrate the relatively high accuracy and reproducibility of cryo-EM models derived by 13 participating teams from four benchmark maps, including three forming a resolution series (1.8 to 3.1 Å). The results permit specific recommendations to be made about validating near-atomic cryo-EM structures both in the context of individual experiments and structure data archives such as the Protein Data Bank. We recommend the adoption of multiple scoring parameters to provide full and objective annotation and assessment of the model, reflective of the observed cryo-EM map density.

qFit 3: Protein and ligand multiconformer modeling for X-ray crystallographic and single-particle cryo-EM density maps.

New X-ray crystallography and cryo-electron microscopy (cryo-EM) approaches yield vast amounts of structural data from dynamic proteins and their complexes. Modeling the full conformational ensemble can provide important biological insights, but identifying and modeling an internally consistent set of alternate conformations remains a formidable challenge. qFit efficiently automates this process by generating a parsimonious multiconformer model. We refactored qFit from a distributed application into software that runs efficiently on a small server, desktop, or laptop. We describe the new qFit 3 software and provide some examples. qFit 3 is open-source under the MIT license, and is available at https://github.com/ExcitedStates/qfit-3.0.

Genomic Predictors of Good Outcome, Recurrence, or Progression in High-Grade T1 Non-Muscle-Invasive Bladder Cancer.

High-grade T1 (HGT1) bladder cancer is the highest risk subtype of non-muscle-invasive bladder cancer with unpredictable outcome and poorly understood risk factors. Here, we examined the association of somatic mutation profiles with nonrecurrent disease (GO, good outcome), recurrence (R), or progression (PD) in a cohort of HGT1 patients. Exome sequencing was performed on 62 HGT1 and 15 matched normal tissue samples. Both tumor only (TO) and paired analyses were performed, focusing on 95 genes known to be mutated in bladder cancer. Somatic mutations, copy-number alterations, mutation load, and mutation signatures were studied. Thirty-three GO, 10 R, 18 PD, and 1 unknown outcome patients were analyzed. Tumor mutational burden (TMB) was similar to muscle-invasive disease and was highest in GO, intermediate in PD, and lowest in R patients (P = 0.017). DNA damage response gene mutations were associated with higher TMB (P < 0.0001) and GO (P = 0.003). ERCC2 and BRCA2 mutations were associated with GO. TP53, ATM, ARID1A, AHR, and SMARCB1 mutations were more frequent in PD. Focal copy-number gain in CCNE1 and CDKN2A deletion was enriched in PD or R (P = 0.047; P = 0.06). APOBEC (46%) and COSMIC5 (34%) signatures were most frequent. APOBEC-A and ERCC2 mutant tumors (COSMIC5) were associated with GO (P = 0.047; P = 0.0002). pT1b microstaging was associated with a genomic cluster (P = 0.05) with focal amplifications of E2F3/SOX4, PVRL4, CCNE1, and TP53 mutations. Findings were validated using external public datasets. These findings require confirmation but suggest that management of HGT1 bladder cancer may be improved via molecular characterization to predict outcome. SIGNIFICANCE: Detailed genetic analyses of HGT1 bladder tumors identify features that correlate with outcome, e.g., high mutational burden, ERCC2 mutations, and high APOBEC-A/ERCC2 mutation signatures were associated with good outcome.

Robust Sequence Determinants of α-Synuclein Toxicity in Yeast Implicate Membrane Binding.

Protein conformations are shaped by cellular environments, but how environmental changes alter the conformational landscapes of specific proteins in vivo remains largely uncharacterized, in part due to the challenge of probing protein structures in living cells. Here, we use deep mutational scanning to investigate how a toxic conformation of α-synuclein, a dynamic protein linked to Parkinson's disease, responds to perturbations of cellular proteostasis. In the context of a course for graduate students in the UCSF Integrative Program in Quantitative Biology, we screened a comprehensive library of α-synuclein missense mutants in yeast cells treated with a variety of small molecules that perturb cellular processes linked to α-synuclein biology and pathobiology. We found that the conformation of α-synuclein previously shown to drive yeast toxicity-an extended, membrane-bound helix-is largely unaffected by these chemical perturbations, underscoring the importance of this conformational state as a driver of cellular toxicity. On the other hand, the chemical perturbations have a significant effect on the ability of mutations to suppress α-synuclein toxicity. Moreover, we find that sequence determinants of α-synuclein toxicity are well described by a simple structural model of the membrane-bound helix. This model predicts that α-synuclein penetrates the membrane to constant depth across its length but that membrane affinity decreases toward the C terminus, which is consistent with orthogonal biophysical measurements. Finally, we discuss how parallelized chemical genetics experiments can provide a robust framework for inquiry-based graduate coursework.

Publisher Correction: The long tail of oncogenic drivers in prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Integrative Molecular Characterization of Resistance to Neoadjuvant Chemoradiation in Rectal Cancer.

PURPOSE: Molecular properties associated with complete response or acquired resistance to concurrent chemotherapy and radiotherapy (CRT) are incompletely characterized.Experimental Design: We performed integrated whole-exome/transcriptome sequencing and immune infiltrate analysis on rectal adenocarcinoma tumors prior to neoadjuvant CRT (pre-CRT) and at time of resection (post-CRT) in 17 patients [8 complete/partial responders, 9 nonresponders (NR)]. RESULTS: CRT was not associated with increased tumor mutational burden or neoantigen load and did not alter the distribution of established somatic tumor mutations in rectal cancer. Concurrent KRAS/TP53 mutations (KP) associated with NR tumors and were enriched for an epithelial-mesenchymal transition transcriptional program. Furthermore, NR was associated with reduced CD4/CD8 T-cell infiltrates and a post-CRT M2 macrophage phenotype. Absence of any local tumor recurrences, KP/NR status predicted worse progression-free survival, suggesting that local immune escape during or after CRT with specific genomic features contributes to distant progression. CONCLUSIONS: Overall, while CRT did not impact genomic profiles, CRT impacted the tumor immune microenvironment, particularly in resistant cases.

Everolimus and pazopanib (E/P) benefit genomically selected patients with metastatic urothelial carcinoma.

BACKGROUND: Metastatic urothelial carcinoma (mUC) is a genomically diverse disease with known alterations in the mTOR pathway and tyrosine kinases including FGFR. We investigated the efficacy and safety of combination treatment with everolimus and pazopanib (E/P) in genomically profiled patients with mUC. METHODS: mUC patients enrolled on a Phase I dose escalation study and an expansion cohort treated with E/P were included. The primary end point was objective response rate (ORR); secondary end points were safety, duration of response (DOR), progression-free survival (PFS) and overall survival (OS). Patients were assessed for mutations and copy number alterations in 300 relevant cancer-associated genes using next-generation sequencing and findings were correlated with outcomes. Time-to-event data were estimated with Kaplan-Meier methods. RESULTS: Of the 23 patients enrolled overall, 19 had mUC. ORR was 21% (one complete response (CR), three partial responses (PR), eight with stable disease (SD). DOR, PFS and OS were 6.5, 3.6, and 9.1 months, respectively. Four patients with clinical benefit (one CR, two PR, one SD) had mutations in TSC1/TSC2 or mTOR and a 5th patient with PR had a FGFR3-TACC3 fusion. CONCLUSIONS: Combination therapy with E/P is safe in mUC and select patients with alterations in mTOR or FGFR pathways derive significant clinical benefit.

Racial disparity in quality of care and overall survival among black vs. white patients with muscle-invasive bladder cancer treated with radical cystectomy: A national cancer database analysis.

OBJECTIVES: To examine the impact of race on quality of care and overall survival (OS) among patients with muscle invasive bladder cancer (MIBC) treated with radical cystectomy (RC) in the U.S. MATERIALS & METHODS: Our cohort consisted of 12,652 patients receiving RC for MIBC within the National Cancer Database from 2004 to 2012. Patients were stratified by race (Black non-Hispanic vs. White non-Hispanic) and imbalances in patient characteristics mitigated using propensity score weighting. Logistic and Cox regressions examined the impact of race on quality of care metrics (receipt of pelvic lymph node dissection (PLND), lymph node count, hospital volume, length of stay, delay of treatment) and on OS. The difference in OS was expressed as Delta, and stratified by facility-type, hospital volume, and region. RESULTS: Blacks were less likely to receive PLND (odds ratio [OR] 0.70, 95% confidence interval [CI]: 0.55-0.91), or to have a greater number of lymph nodes removed (OR 0.76, 95%CI: 0.64-0.90). They exhibited greater length of stay (OR 1.34, 95%CI: 1.13-1.59), and delay of RC among recipients of neoadjuvant chemotherapy (OR 2.59, 95%CI: 1.77-3.85) (all P ≤ 0.001). Notably, utilization of neoadjuvant chemotherapy in advanced disease stages was more common in blacks (OR 2.82, 95%CI: 1.93-4.13, P < 0.001). Additionally, Black race was associated with inferior OS (Hazard ratio 0.87, 95%CI: 0.79-0.97, P < 0.014). Disparities in OS varied based on facility type and geographical region, but not hospital volume. Specifically, Blacks had worse OS when treated in a community cancer program (Delta 0.42, 95%CI: 0.28-0.57,P < 0.001), or within New England/Middle Atlantic region (Delta 0.16, 95% CI: 0.07-0.24,P < 0.001). CONCLUSION: Black race is an independent predictor of inferior quality of care and OS in patients undergoing RC for MIBC. Survival disparities vary based on geographical region and facility type. Notably, the OS disparity appears to have narrowed in comparison to previous studies.

Structural Alterations Driving Castration-Resistant Prostate Cancer Revealed by Linked-Read Genome Sequencing.

Nearly all prostate cancer deaths are from metastatic castration-resistant prostate cancer (mCRPC), but there have been few whole-genome sequencing (WGS) studies of this disease state. We performed linked-read WGS on 23 mCRPC biopsy specimens and analyzed cell-free DNA sequencing data from 86 patients with mCRPC. In addition to frequent rearrangements affecting known prostate cancer genes, we observed complex rearrangements of the AR locus in most cases. Unexpectedly, these rearrangements include highly recurrent tandem duplications involving an upstream enhancer of AR in 70%-87% of cases compared with <2% of primary prostate cancers. A subset of cases displayed AR or MYC enhancer duplication in the context of a genome-wide tandem duplicator phenotype associated with CDK12 inactivation. Our findings highlight the complex genomic structure of mCRPC, nominate alterations that may inform prostate cancer treatment, and suggest that additional recurrent events in the non-coding mCRPC genome remain to be discovered.

The Clinical Activity of PD-1/PD-L1 Inhibitors in Metastatic Non-Clear Cell Renal Cell Carcinoma.

Programmed death 1 (PD-1) and PD ligand 1 (PD-L1) inhibitors have shown activity in metastatic clear cell renal cell carcinoma (ccRCC). Data on the activity of these agents in patients with non-clear cell RCC (nccRCC) or patients with sarcomatoid/rhabdoid differentiation are limited. In this multicenter analysis, we explored the efficacy of PD-1/PD-L1 inhibitors in patients with nccRCC or sarcomatoid/rhabdoid differentiation. Baseline and follow-up demographic, clinical, treatment, and radiographic data were collected. The primary endpoint was objective response rate. Secondary endpoints include time-to-treatment failure (TTF), overall survival (OS), and biomarker correlates. Forty-three patients were included: papillary (n = 14; 33%), chromophobe (n = 10; 23%), unclassified (n = 9; 21%), translocation (n = 3; 7%), and ccRCC with sarcomatoid differentiation (n = 7, 16%). Of those 43 patients, 11 patients (26%) had sarcomatoid and/or rhabdoid differentiation (n = 7 with ccRCC; n = 4 nccRCC). Overall, 8 patients (19%) objectively responded, including 4 patients (13%) who received PD-1/PD-L1 monotherapy. Responses were observed in patients with ccRCC with sarcomatoid and/or rhabdoid differentiation (n = 3/7, 43%), translocation RCC (n = 1/3, 33%), and papillary RCC (n = 4/14, 29%). The median TTF was 4.0 months [95% confidence interval (CI), 2.8-5.5] and median OS was 12.9 months (95% CI, 7.4-not reached). No specific genomic alteration was associated with clinical benefit. Modest antitumor activity for PD-1/PD-L1-blocking agents was observed in some patients with nccRCC. Further prospective studies are warranted to investigate the efficacy of PD-1/PD-L1 blockade in this heterogeneous patient population. Cancer Immunol Res; 6(7); 758-65. ©2018 AACR.

Targeted genomic landscape of metastases compared to primary tumours in clear cell metastatic renal cell carcinoma.

BACKGROUND: The genomic landscape of primary clear cell renal cell carcinoma (ccRCC) has been well described. However, little is known about cohort genomic alterations (GA) landscape in ccRCC metastases, or how it compares to primary tumours in aggregate. The genomic landscape of metastases may have biological, clinical, and therapeutic implications. METHODS: We collected targeted next-generation sequencing mutation calls from two independent cohorts and described the metastases GA landscape and descriptively compared it to the GA landscape in primary tumours. RESULTS: The cohort 1 (n = 578) consisted of 349 primary tumours and 229 metastases. Overall, the most common mutations in the metastases were VHL (66.8%), PBRM1 (41.87%), and SETD2 (24.7%). TP53 was more frequently mutated in metastases compared to primary tumours (14.85% versus 8.9%; p = 0.031). No other gene had significant difference in the cohort frequency of mutations between the metastases and primary tumours. Mutation burden was not significantly different between the metastases and primary tumours or between metastatic sites. The second cohort (n = 257) consisted of 177 primary tumours and 80 metastases. No differences in frequency of mutations or mutational burden were observed between primaries and metastases. CONCLUSIONS: These data support the theory that ccRCC primary tumours and metastases encompass a uniform distribution of common genomic alterations tested by next-generation sequencing targeted panels. This study does not address variability between matched primary tumours and metastases or the change in genomic alterations over time and after sequential systemic therapies.

The long tail of oncogenic drivers in prostate cancer.

Comprehensive genomic characterization of prostate cancer has identified recurrent alterations in genes involved in androgen signaling, DNA repair, and PI3K signaling, among others. However, larger and uniform genomic analysis may identify additional recurrently mutated genes at lower frequencies. Here we aggregate and uniformly analyze exome sequencing data from 1,013 prostate cancers. We identify and validate a new class of E26 transformation-specific (ETS)-fusion-negative tumors defined by mutations in epigenetic regulators, as well as alterations in pathways not previously implicated in prostate cancer, such as the spliceosome pathway. We find that the incidence of significantly mutated genes (SMGs) follows a long-tail distribution, with many genes mutated in less than 3% of cases. We identify a total of 97 SMGs, including 70 not previously implicated in prostate cancer, such as the ubiquitin ligase CUL3 and the transcription factor SPEN. Finally, comparing primary and metastatic prostate cancer identifies a set of genomic markers that may inform risk stratification.

Upper Tract Urothelial Carcinomas: Prognostic Factors and Outcomes in Patients With Non-Lymph Node Distant Metastasis.

BACKGROUND: Upper tract urothelial carcinomas (UTUCs) are increasingly recognized as separate malignancies. Additional insight into clinical outcomes and key prognostic factors are needed. OBJECTIVES: To detail outcomes of patients with UTUCs recurring after radical nephroureterectomy (RNU) and to determine a risk score that predicts outcomes of patients with non-lymph node distant metastasis. DESIGN, SETTING, AND PARTICIPANTS: Chart review of all patients who had an extraurothelial recurrence after RNU for UTUC at Dana-Farber Cancer Institute between 2009 and 2014. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Median overall survival defined as time from chemotherapy for distant relapse to death. Prognostic relevance of performance status, hemoglobin, and receipt of cisplatin were assessed by Cox regression model. RESULTS AND LIMITATIONS: A total of 102 patients were identified, 57 of whom had non-lymph node distant metastases at relapse; 45 received chemotherapy. Median follow-up was 29.8 months; median overall survival was 14.7 months. Objective response rate to any chemotherapy in the first-line setting was only 22%. Hemoglobin > 11 g/dL and receipt of cisplatin was associated with numerically longer median survival but did not reach statistical significance in univariate and multivariate analysis. Prognostic risk score scale including hemoglobin < 11 g/dL and receipt of cisplatin was inversely associated with survival, with scores of 0, 1, and 2 leading to median survival of 19.0, 14.9, and 7.2 months (P = .38), respectively. CONCLUSIONS: Advanced UTUC portends a poor prognosis, and most patients cannot receive cisplatin-based chemotherapy. A risk score that includes anemia and receipt of cisplatin helps stratify patients with distant metastasis for inclusion into eventual clinical trials. More studies are needed to validate these findings. PATIENT SUMMARY: Metastatic UTUC is an aggressive disease, where anemia and ineligibility to receive cisplatin are adverse features associated with shorter survival.

Neoadjuvant chemotherapy prior to radical cystectomy for muscle-invasive bladder cancer with variant histology.

BACKGROUND: Neoadjuvant chemotherapy in pure urothelial bladder cancer provides a significant survival benefit. However, to the authors' knowledge, it is unknown whether this benefit persists in histological variants. The objective of the current study was to assess the effect of neoadjuvant chemotherapy on the probability of non-organ-confined disease and overall survival after radical cystectomy (RC) in patients with histological variants. METHODS: Querying the National Cancer Data Base, the authors identified 2018 patients with histological variants who were undergoing RC for bladder cancer between 2003 and 2012. Variants were categorized as micropapillary or sarcomatoid differentiation, squamous cell carcinoma, adenocarcinoma, neuroendocrine tumors, and other histology. Logistic regression models estimated the odds of non-organ-confined disease at the time of RC for each histological variant, stratified by the receipt of neoadjuvant chemotherapy. Cox regression models were used to examine the effect of neoadjuvant chemotherapy on overall mortality in each variant subgroup. RESULTS: Patients with neuroendocrine tumors (odds ratio [OR], 0.16; 95% confidence interval [95% CI], 0.08-0.32 [P<.001]), micropapillary differentiation (OR, 0.30; 95% CI, 0.10-0.95 [P=.041]), sarcomatoid urothelial carcinoma (OR, 0.40; 95% CI, 0.17-0.94 [P=.035]), and adenocarcinoma (OR, 0.24; 95% CI, 0.06-0.91 [P=.035]) were less likely to harbor non-organ-confined disease at the time of RC when treated with neoadjuvant chemotherapy. An overall survival benefit for neoadjuvant chemotherapy was only found in patients with neuroendocrine tumors (hazard ratio, 0.49; 95% CI, 0.33-0.74 [P=.001]). CONCLUSIONS: Patients with neuroendocrine tumors benefit from neoadjuvant chemotherapy, as evidenced by better overall survival and lower rates of non-organ-confined disease at the time of RC. For tumors with micropapillary differentiation, sarcomatoid differentiation, or adenocarcinoma, neoadjuvant chemotherapy decreased the frequency of non-organ-confined disease at the time of RC. However, this favorable effect did not translate into a statistically significant overall survival benefit for these patients, potentially due to the aggressive tumor biology. Cancer 2017;123:4346-55. © 2017 American Cancer Society.

Differential Expression of PD-L1 in High Grade T1 vs Muscle Invasive Bladder Carcinoma and its Prognostic Implications.

PURPOSE: PD-L1 is expressed on tumor cells and tumor immune cell infiltrates. In metastatic bladder cancer increased tumor immune cell infiltrate PD-L1 positivity correlated with better overall survival. However, to our knowledge in high grade T1 bladder tumors positivity on tumor cells and tumor immune cell infiltrates, and correlation with outcomes or pathological features remain unknown. MATERIALS AND METHODS: Formalin fixed, paraffin embedded tumor samples from 140 patients with clinically annotated, high grade T1 bladder tumors were retrieved. All patients were initially diagnosed with high grade T1 bladder tumors by transurethral resection, subsequently received bacillus Calmette-Guérin and had a median followup of 7.4 years. PD-L1 positivity on initial transurethral resection was evaluated by immunohistochemistry using a mouse monoclonal antiPD-L1 antibody (405.9A11). Tumor cell PD-L1 positivity was defined as staining of 5% of the tumor cell membrane. Tumor immune cell infiltrate PD-L1 positivity was scored based on the extent of infiltrate and the percent of positive cells. The Fisher exact test was used to assess associations of PD-L1 positivity with disease outcomes, carcinoma in situ presence and the difference between high grade T1 bladder tumors and muscle invasive bladder cancer. RESULTS: Among 140 patients with high grade T1 bladder tumors tumor cells and tumor immune cell infiltrate PD-L1 positivity was seen in 6 (4%) and 48 (34.3%), respectively. In a subset of 106 patients with adequate followup PD-L1 positivity did not correlate with disease outcomes on tumor cells (p = 0.3) or on tumor immune cell infiltrates (p = 0.47). PD-L1 positivity also did not correlate with the presence of carcinoma in situ. Tumor cell PD-L1 positivity was significantly less in high grade T1 bladder tumors than in muscle invasive bladder cancer (p <0.001). CONCLUSIONS: PD-L1 is widely expressed on tumor immune cell infiltrates but not on tumor cells in high grade T1 bladder tumors. We did not find a correlation between PD-L1 positivity and outcomes or carcinoma in situ presence. Tumor cell PD-L1 positivity is significantly lower in high grade T1 bladder tumors than in muscle invasive bladder cancer.